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Abstract 

There is a growing concern for the health and safety of commercial aircrew and 
passengers due to their exposure to ionizing radiation with high linear energy transfer 
(LET), particularly at high latitudes. The International Commission of Radiobiological 
Protection (ICRP), the EPA, and the FAA consider the crews of commercial aircraft as 
radiation workers. During solar energetic particle (SEP) events, radiation exposure can 
exceed annual limits, and the number of serious health effects is expected to be quite high 
if precautions are not taken. There is a need for a capability to monitor the real-time, 
global background radiations levels, from galactic cosmic rays (GCR), at commercial 
airline altitudes and to provide analytical input for airline operations decisions for altering 
flight paths and altitudes for the mitigation and reduction of radiation exposure levels 
during a SEP event. The Nowcast of Atmospheric Ionizing Radiation for Aviation Safety 
(NAIRAS) model is new initiative to provide a global, real-time radiation dosimetry 
package for archiving and assessing the biologically harmful radiation exposure levels at 
commercial airline altitudes. The NAIRAS model brings to bear the best available suite 
of Sun-Earth observations and models for simulating the atmospheric ionizing radiation 
environment. Observations are utilized from ground (neutron monitors), from the 
atmosphere (the METO analysis), and from space (NASA/ACE and NOAA/GOES). 
Atmospheric observations provide the overhead shielding information and the ground- 
and space-based observations provide boundary conditions on the GCR and SEP energy 
flux distributions for transport and dosimetry simulations. Dose rates are calculated using 
the parametric AIR (Atmospheric Ionizing Radiation) model and the physics-based 
HZETRN (High Charge and Energy Transport) code. Empirical models of the near-Earth 
radiation environment (GCR/SEP energy flux distributions and geomagnetic cut-off 
rigidity) are benchmarked against the physics-based CMIT (Coupled Magnetosphere- 
Ionosphere-Thermosphere) and SEP-trajectory models. 

1.0 Introduction 
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1. 1 Scientific Background 

Atmospheric ionizing radiation is of interest to air transportation safety assessment 
because it’s the primary source of human exposure to radiations with high linear energy 
transfer (LET). High-LET radiation is effective at producing chemically active radicals in 
biological tissues that alter the cell function or result in cell death. Consequently, there is 
increased concern for potential health outcomes among passengers and crew in 
commercial aviation [ Wilson et al., 2003]. Atmospheric ionizing radiation is produced by 
extraterrestrial radiations incident on the Earth’s atmosphere. There are two sources of 
the extraterrestrial radiations: the ever-present, background galactic cosmic rays (GCR), 
with origins outside the solar system, and transient solar energetic particles (SEP) 
associated with solar storm activity lasting several hours to days with widely varying 
intensity. 

GCR consist of roughly 90% protons and 8% helium nuclei with the remainder being 
heavier nuclei and electrons [Gaisser, 1990]. When these particles penetrate the magnetic 
fields of the solar system and the Earth and reach the Earth’s atmosphere, they collide 
with air molecules and create cascades of secondary radiations of every kind [Reitz, 

1993]. The collisions are primarily due to Coulomb interactions of the GCR particle with 
orbital electrons of the air molecules, delivering small amounts of energy to the orbital 
electrons and leaving behind electron-ion pairs [Wilson et al., 1991]. The ejected 
electrons usually have sufficient energy to undergo similar ionizing events. The cosmic 
ions lose a small fraction of their energy and must suffer many collisions before slowing 
down. On rare occasions the cosmic ion will collide with the nucleus of an air atom in 
which large energies are exchanged and the ion and nucleus are dramatically changed by 
the violence of the event. The remnant nucleus is highly disfigured and unstable, emitting 
further air nuclear constituents and decaying through the usual radioactivity channels 
[Wilson et al, 1991]. The most important secondary particle created in GCR-air 
interactions is the neutron. Because of its charge neutrality, the neutron penetrates deep 
into the atmosphere, causing further ionization events along its path and contributing over 
half the atmospheric radiation exposure [Wilson et al., 1993]. Furthermore, neutron 
exposures pose a relatively high health risk, since the massive low-energy ions resulting 
from neutron interactions always produce copious ions in the struck cell and repair is less 
efficient for these events [Wilson et al., 2000]. 

The intensity of the atmospheric radiations, composed of GCR primary and secondary 
particles, their energy distribution, and their effects on aircraft occupants vary with 
altitude, location in the geomagnetic field, and the time in the sun’s magnetic activity 
(solar) cycle [Reitz, 1993; Wilson, 2000; and Heinrich et al., 1999]. The atmosphere 
provides shielding, which depends on the overhead atmospheric depth. The geomagnetic 
field provides a different kind of shielding, by deflecting low-momentum charged 
particles back to space. Because of the orientation of the geomagnetic field, which is 
predominately dipolar in nature, the polar regions are susceptible to penetrating GCR 
(and SEP) particles. At each geographic location, the minimum momentum per unit 
charge (magnetic rigidity) a vertically incident particle can have and still reach a given 
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location above the earth is called the geomagnetic vertical cutoff rigidity. The local flux 
of incident GCR at a given time varies widely with geomagnetic location and the solar 
modulation level. When solar activity is high, GCR flux is low, and vice versa. The 
dynamical balance between outward convective flux of solar wind and the inward 
diffusive flux of GCR is responsible for the anti-correlation between the incident GCR 
and the level of solar activity [Clem et al.; 1996; Parker, 1965]. 

It is now generally understood that SEP events arise from coronal mass ejections (CME) 
from active regions of the solar surface [Kahler, 2001; Wilson et al., 2004]. The CME 
propagates through interplanetary space carrying along with it the local surface magnetic 
field frozen into the ejected mass. There is a transition (shock) region between the nonnal 
sectored magnetic structure of interplanetary space and the fields frozen into the ejected 
mass, which forms a transition region (shock) where the interplanetary gas is accelerated 
forming the SEP. As the accelerated region passes an observation point, the flux intensity 
is observed to increase dramatically, and no upper limit in intensity is known within the 
shock region. The SEP energy spectrum obtained in the acceleration process is related to 
the plasma density and CME velocity. During a solar storm CME event, the number flux 
and energy flux distributions incident at Earth’s atmosphere are a combination of the 
GCR and SEP distributions. The SPE-air interaction mechanisms are the same as GCR- 
air interactions described above. The atmospheric radiations caused by a SEP also vary 
with altitude and geomagnetic field. 

1.2 Identification of Management/Policy Issues and End-User Community 

GCR radiations that penetrate the atmosphere and reach the ground are low intensity. 
However, the intensities are more than two orders of magnitude greater at commercial 
aircraft altitudes. At the higher altitudes of High Speed Civil Transport (HSCT), the GCR 
intensity is another two orders of magnitude higher [ Wilson et al., 2003], When the 
possibility of high-altitude supersonic commercial aviation was first seriously proposed 
(The Supersonic Transport program proposed in 1961), Foelsche brought to light a 
number of concerns about associated atmospheric radiation exposure due to GCR and 
SEP, including the secondary radiations [Foelsche, 1961; Foelsche and Graul, 1962], 
Subsequently, Foelsche et al. [1974] conducted a detailed study of atmospheric ionizing 
radiation at high altitudes from 1965 to 1971 at the NASA Langley Research Center 
(LaRC). The study included a comprehensive flight program in addition to theoretical 
investigations. The measured data and theoretical calculations were integrated into a 
parametric Atmospheric Ionizing Radiation (AIR) model [Wilson et al, 1991]. Prior to 
that study the role of atmospheric neutrons in radiation exposure was generally regarded 
as negligible [Upton et al., 1966]. The LaRC studies revealed neutron radiation to be a 
major contributor to aircraft GCR exposure. Still the exposure levels were comfortably 
below allowable exposure limits for the block hours typical of airline crews of that time, 
except during a possible SEP event (less than 500 block hours were typical of the 1960’s, 
although regulation allowed up to 1000 hours). 

There have been a number of significant changes since the original work of Foelsche 
[Wilson, 2002], A partial list of these changes, relevant to the development of the 
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NAIRAS model, are: (1) the highly ionizing components of atmospheric radiations are 
found to be more biologically damaging than previously assumed and the associated 
relative biological effectiveness for fatal cancer has been increased [ICRU 1986; ICRP 
1991]; (2) recent studies on developmental injury in mice embryos indicate large relative 
biological effectiveness for protection in prenatal exposures [Jiang et a/., 1994]; (3) 
recent epidemiological studies (especially the data on solid tumors) and more recent 
atom-bomb survivor dosimetry have resulted in higher radiation risk coefficients for 
gamma rays [UNSCEAR 1988; NAS/NRC 1980; ICRP 1991], resulting in lower proposed 
permissible limits [ ICRP 1991; NCRP 1993]; (4) subsequent to deregulation of the airline 
industry, flight crews are logging greatly increased hours [Bramlitt, 1985; Wilson and 
Townsend, 1988; Friedberg et al., 1989; Banish, 1990]; and (5) airline crew members are 
now classified as radiation workers [McMeekin, 1990; ICRP 1991]. 

The last point (i.e., (5)) is particularly illuminating. The International Commission on 
Radiological Protection (ICRP), as well as the EPA and FAA, consider the crews of 
commercial aircraft as radiation workers [ Wilson et al., 2003]. The FAA estimates 
annual subsonic aircrew exposures to range from 0.2 to 9.1 mSv compared to 0.5 mSv 
exposure of the average nuclear power plant worker. Aircrews may receive exposures 
above recently recommended allowable limits for even radiation workers when flying the 
maximum allowable number of flight hours. Although as a group the health risks of 
aircrew are low, Band [1990] found increased risks of several types of cancer among 
Canadian commercial pilots. There is further concern for prenatal injury in high altitude 
flight, as the US National Institute for Occupational Safety and Health continues to study 
early pregnancy outcomes among commercial flight attendants [Grajewski et al., 1994; 
Whelan, 2002]. Frequent- flyer business passengers are likely exposed to even higher 
doses than aircrew, since flight hours are not restricted for airline passengers. In addition, 
if a large SEP occurs during flight, both passengers and crew may greatly exceed 
allowable limits and potentially serious health outcomes are possible [ Banish , 2004]. 

There are currently no data or models that can effectively map the atmospheric radiation 
field on a continuous basis in order to track and regulate radiation exposure and 
associated health risks to aircrew and passengers on commercial flights. 

1.3 Model Objectives, Benefits, and Relevance 

Recognizing the potential impact on present day passenger and crew exposures, due to 
the changes since the original work of Foelsche, as described above, further studies were 
started at LaRC. The resulting flight package was a collaboration of fourteen institutions 
in five countries and consisted of eighteen instruments. New measurements [ Wilson et 
al., 2003] were made and new advances in theoretical modeling [ Clem et al., 1996], 
which culminated in the AIR workshop [ Wilson et al., 2003]. 

Following the recent LaRC-sponsored AIR workshop, a number of recommendations for 
future work were put forth [ Wilson et al., 2003]. The recommendations relevant to the 
NAIRAS model development are: (1) utilize satellite input data to provide real-time 
mapping of GCR and SEP radiation levels to provide guidance in exposure avoidance; 
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and (2) utilize state-of-the-art transport codes and nuclear databases to generate input data 
to the AIR model. The need to combine satellite observations of the meteorological fields 
and the space environment variables with HZE (high charge and energy) particle 
transport codes was further highlighted in a recent Airline Space Weather Workshop 
[ Friedberg , 2004], The objective of the NAIRAS model is to address these two 
recommendations, but significantly go beyond recommendation (2) by using physics- 
based, state-of-the-art transport code directly in simulating the atmospheric radiation 
exposure levels. 

NAIRAS is a new initiative to develop a prototype, global, Nowcast of Atmospheric 
Ionizing Radiation for Aviation Safety (NAIRAS) model for calculating dose equivalent 
rates at commercial airline altitudes. The NAIRAS results will provide a continuous 
assessment of the atmospheric ionizing radiation field needed for the commercial airlines 
to track individual aircrew radiation exposures levels, in order that the airlines and the 
FAA can develop policy and procedures for radiation exposure limits and exposure 
mitigation to aircrew. During SEP events, NAIRAS will provide timely data output 
necessary for airline management to make critical decisions that balance cost to flight 
path alterations against radiation exposure and health risks to passenger and crew. 

The most significant and innovative features of the NAIRAS model are: (1) the 
atmospheric transport is simulated using HZETRN (High Charge and Energy Transport), 
a state-of-the-art, physics-based HZE+neutron+meson+muon transport code, which is (2) 
driven by real-time measurements of the solar-terrestrial environment - i.e., 
meteorological data of atmospheric density, observation-based models of GCR/SEP 
differential number flux spectrum (DNFS), and observation-based models of 
geomagnetic cutoff rigidity. Currently NOAA/SEC maintains a web site with space 
weather forecast products for the aviation community. The products include 
measurements of several space radiation environment parameters and a forecast on the 
likelihood and the expected level of space weather activity. The NOAA aviation products 
do not include the radiation fields that effect human health. The FAA’s Civil Aerospace 
Medical Institute provides a web interface to the CARI-6 program. CARI-6 calculates 
effective dose rates of the GCR radiations for user-defined flight paths. The FAA web 
site does not maintain a real-time, global database. 

There are no existing data or models that provide a comprehensive (i.e., comprehensive 
in terms of input observation data included and comprehensive in terms of the transport 
physics included in the real-time calculations), global, real-time assessment of the 
radiation fields that affect human health and safety. Thus, the NAIRAS model concept 
provides an atmospheric radiation exposure assessment that significantly extends current 
capabilities. 

2.0 Description of the Parametric AIR Model 

The original LaRC study (1965 to 1971) commissioned over 300 flights over most of the 
duration of solar cycle 20 on high-altitude aircraft and balloons to study both the 
background radiation levels over the solar cycle and to make measurements during SEP 
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events. The LaRC flight package consisted of a 1-10 MeV neutron spectrometer, tissue 
equivalent ion chamber, and nuclear emulsion for nuclear reaction rates in tissue. Monte 
Carlo calculations [ Wilson et al, 1970; Lambiotte et al., 1971] for incident GCR protons 
were used to extend the neutron spectrum to high energies. The measured data was 
combined with the theoretical calculations and integrated into a parametric Atmospheric 
Ionizing Radiation (AIR) model, parameterized by neutron monitor count rate, 
geomagnetic vertical cutoff rigidity, and atmospheric depth. Solar cycle modulation of 
the GCR spectrum is parameterized by the ground-level neutron monitor count rates. 
Geomagnetic momentum shielding and overhead atmospheric shielding are 
parameterized by the vertical geomagnetic cutoff rigidity and atmospheric depth, 
respectively. The neutron flux (cm'" sec' ) component to the atmospheric radiations is 
converted to dose equivalent and total dose using 3.14 pSv cm" sec hr' and 0.5 pGy cm' 
sec hr' 1 , respectively. The charged particle component to the atmospheric radiations is 
obtained from data taken by Neher [1961. 1967, 1971] and Neher and Anderson [1962] 
as compiled S. B. Curtis (Boeing 1969) and utilized by Wallace and Sondhaus [1978]. 

The charge particle atmospheric ionization rates are directly converted to dose equivalent 
and total dose using measurement data from the tissue equivalent ion chamber. Nuclear 
stars in tissue are estimated from the nuclear emulsion measurement data after subtraction 
of the neutron-induced stars [ Wilson et al., 1991]. 

Considerable progress was made in the recent LaRC study [ Wilson et al., 2005a, 2003a], 
In particular, improvements were made in the high-energy neutron spectrum from a 
combination of flight measurements and new theoretical calculations using the FLUKA 
transport code [ Clem et al., 2003]. Furthermore, an improved model was developed for 
geomagnetic cutoff rigidity applicable for years 1945 to 2020 which can also incorporate 
geomagnetic storm effects [De Ange/is et al., 2003], 

Figure 1 shows altitude profiles of neutron flux and ionization rates computed from the 
AIR model for summer and winter atmospheric conditions at various latitudes in the 
northern hemisphere, for both solar maximum and solar minimum conditions. Figure 2 
shows the corresponding profiles of dose equivalent rates. The two most noticeable 
features are: (1) the significant increase in flux, ionization rates, and dose-equivalent rates 
at high-latitudes, and (2) the peak in these quantities occur near the typical cruising 
altitudes of commercial aircraft flying international routes (~ 10-12 km). The low altitude 
results are less reliable because of the limited altitude range of the balloon and flight 
measurements used to develop the AIR parameterization. 

Figures 1 and 2 show the dose-equivalent rates for summer and winter northern 
hemispheres at 12 km for both solar minimum and solar maximum conditions, 
respectively. The North Atlantic flight corridor is one of the busiest in the world and it is 
among the most highly exposed routes in airline operations. Flights over Canada are 
among the most highly exposed. Much of European flight is subject to somewhat lower 
exposure levels. From Table 1 it is clear that aircrew flying the Northern Atlantic or 
Canadian routes can exceed allowable annual exposure levels (10 mSv, see footnote b in 
Table 1) in a 1000-hour block at solar minimum. For solar maximum condition, aircrew 
can reach 60% to 70% of the annual recommend allowance in a 1000-hour block. The 
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occurrence of a SEP could increase radiation exposure well over recommended and 
allowable levels. 

The AIR model will be used as an intermediate tool to develop the interface between 
real-time neutron monitor data, the atmospheric depth data, and the observation-based 
geomagnetic cutoff model. This will allow simultaneous development of the integration 
of the components of the HZETRN code and the interface between HZETRN and the 
data-driven GCR/SEP models. Once the geomagnetic cutoff model and the atmospheric 
depth data has been validated, verified, and benchmarked using the AIR model, these 
models and data input will be to integrated into HZETRN. A common I/O interface and 
data definitions will make this step effortless. Thus, HZETRN will replace the AIR 
model for dose rate calculations, since the physics-based HZETRN is the best alternative. 
The purpose of the AIR model is to facilitate the integration, testing, and benchmarking 
of the components of the NAIRAS model. 

3.0 Description of the NAIRAS Model Components 

3. 1 Transport Code 

The LaRC HZETRN transport code has a long and successful history for rapidly and 
accurately modeling the particle radiation fields in the space environment. The LaRC 
code is used to calculate dosimetry parameters on the International Space Station (ISS) 
and assess astronaut risk to space radiations, including spacecraft and human geometry 
for final exposure evaluation. It is used to develop design tools for materials research for 
radiation shielding protection, to calculate HZE propagation through Earth’s atmosphere, 
and to evaluate radiation exposures for epidemiological studies [ Wilson et al., 1997, 

2003, and references therein]. 

The relevant transport equations are the linear Boltzmann equations derived on the basis 
of conservation principles [ Wilson et al., 1991]. The transport equation for the flux 
density (fj(x, Q. E) for particle type j is given by 


o • v <t>j (x, n, E) = X J o> E')<t> k (X, SI', E') - cj j (SI, E)(f> j (x,Sl,E) ( 1 ) 

k 

where o](E) and <Jjk(Q,Q',E,E) are the target medium macroscopic cross sections. The 
<jjk(Q£2',E,E) represent all those processes by which type k particles moving in direction 
Q' with energy E' produce a type j particle in direction Z? with energy E (including decay 
processes). The total cross section <jj(Sl,E) with the target medium for each particle 

type is 


a J (n,E) = a Mt (Sl,E) + cr jel (Sl,E) + a jr (Q, E). (2) 

The first term above refers to collisions with atomic electrons, the second tenn refers to 
elastic ion-nucleus scattering, and the third contains all relevant nuclear reactions and 
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radioactive decay processes. The corresponding differential cross sections are similarly 
ordered. 

The solution of (1) involves hundreds of multidimensional integral-differential equations, 
which are coupled together by thousands of cross terms and must be solved self- 
consistently subject to boundary conditions ultimately related to the external 
environment. HZETRN determines the solution of (1) for the transport of HZE (High 
Charge and Energy) particles, light-ions (protons through alpha particles), neutrons, and 
mesons and muons. The transport of each particle type described above employs different 
approximations and requires different solution approaches. Details of the analytical and 
computational solution approaches implemented in HZETRN are given by Wilson et al. 
[2004b, 1997, 1991, and references therein], C lowdsley et al [2000, 2002], and Blattnig 
et al. [2004, 2005]. 

Figures 5-8 show examples of atmospheric GCR and SEP flux and dose-equivalent rates 
computed by HZETRN for the solar minimum conditions and the September 1989 SEP 
event. Figure 5 shows the differential flux for neutrons and selected ions for solar 
minimum at various atmospheric depths. The neutron flux quickly builds up and even 
exceeds the proton flux. At depths greater than 100 g/cm2, the heavy ion flux quickly 
approaches zero, leaving mostly protons, neutrons, electrons and muons. Figure 6 shows 
the corresponding total dose-equivalent rates at various atmospheric depths for solar 
minimum. The event-integrated proton and neutron flux and total dose-equivalent rates 
for the September 1989 SEP events are shown in Figures 7 and 8, respectively. 

The next sections briefly discuss the models that specify the input data for the HZETRN 
transport calculations, namely, the GCR/SEP flux incident on Earth’s atmosphere, the 
geomagnetic cutoff rigidity, and atmospheric pressure versus altitude. The results 
presented in Figures 5-8 assumed a zero cutoff rigidity for all incident ions. 

3.2 GCR Model 

The solutions to the Boltzmann transport equation (1) are unique in any convex region for 
which the inbound flux of each particle type is specified everywhere on the bounding 
surface [ Wilson et al., 1997]. For real-time transport calculations of GCR particles, we 
use the Badhwar and O ’Neill [1996, 1994, 1993, 1992, 1991] model to specify the 
incident GCR DNFS at the top of the atmosphere. The GCR DNFS is derived by solving 
a steady-state Fokker-Planck equation for the transport of GCR particles through the 
heliosphere. The local interstellar spectrum (LIS) at 100 AU is parameterized by 

j LIS (E) = jJ s (E + E o r. (3) 

In the above equation, E and E 0 are the particle’s kinetic and rest energy per nucleon, 
respectively, and /? is the particle’s velocity divided by the speed of light. The free 
parameters ( j Q , S , and y ) are detennined by a fit to the NASA/ ACE data, as described 
below. 
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The diffusion coefficient in the Fokker-Planck transport equations are parameterized by 

k(r, t ) = (k o / V sw )J3R[1 + (r / r ) 2 ] / 0(t). (4) 

where V sw is the solar wind velocity, /? is the particle’s velocity divided by the speed of 
light, R is the particle’s rigidity, radial distances from the sun in units of AU are r and r 0 , 
® is the heliospheric potential, and t is time. Nominal values are specified for k 0 , V sw , 
and r 0 . The single fit parameter is taken to be the heliospheric potential. 

Figure 9 shows a sample output for the Badhwar and O’Neill GCR model for a select 
number of ions. The figure shows the LIS spectra and the spectra at solar maximum and 
solar minimum conditions. The proton and alpha spectra were fit to IMP-8 data. The 
lithium through nickel (Z=3-28) spectra were fit to CRIS measurements taken from 
NASA/ACE satellite (50-500 MeV/nucleon). The high-energy spectra (1-35 
GeV/nucleon) were fit to C2 observations measured by the NASA/HEAO-3 satellite 
[Engelmann et al., 1999]. 

The recently updated version of the Badhwar and O’Neill GCR model can be driven by 
ground-based neutron monitor count rate measurements. A reference heliospheric 
potential was fit to the ACE/CRIS oxygen spectra. The free parameters in (3) were fit 
using the reference heliospheric potential and ACE/CRIS spectral measurements for the 
remaining nuclei. A simple linear scaling relationship was determined between the 
heliospheric potential derived from the ACE/CRIS measurements and the CLIMAX 
neutron monitor count rate, which enables the heliospheric potential to be calculated from 
1951 to the present. Thus, the real-time modulation of the GCR DNFS by the solar wind 
can be captured using neutron monitor data. Real-time neutron monitor data, with 
geomagnetic cutoff rigidities of 2 GV or greater, will be obtained from sites at IZMIRAN 
(Moscow, Russia), YAKUTSK (Russia), and LOMNICKY (Slovakia). We will derive 
linear scaling formulas between the real-time neutron monitor data and the heliospheric 
potential derived during the ACE/CRIS measurement period. 

3.3 SEP Spectrum 

SEP proton and alpha DNFS will be obtained in real-time from NOAA’s Geostationary 
Operational Environmental Satellite (GOES) Space Environment Monitor (SEC) 
measurements. The Energetic Particle Sensor (EPS) and the High Energy Proton and 
Alpha Detector (HEP AD) sensors on GOES/SEC measure energetic differential proton 
and alpha flux. EPS provides seven-channel differential proton flux from 0.8 to 500 MeV 
and six-channel differential alpha flux from 4 to 500 MeV per nucleon. HEPAD extends 
the EPS energy ranges to greater than 700 MeV for protons and 3400 MeV per nucleon 
for alpha particles. These measurements will be used to fit a power law spectrum [Wilson 
et al., 2003a] for the incident proton/alpha flux such that the DNFS integrated over the 
spectral channel subintervals agree with the EPS and HEPAD measurements, in the least- 
square sense. 
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3.4 Geomagnetic Cutoff Rigidity 

The geomagnetic cutoff rigidity detennines the minimum energy for the transport of 
GCR/SEP particles through the atmosphere using HZETRN. A baseline semi-empirical 
cutoff rigidity model is based on the worldwide grid maps produced by Shea and Smart 
[1983] and Smart and Shea [2000, 2001]. Vertical cutoff rigidities from Shea and Smart 
data are shown in Figure 10. Note the correlation between high rigidity in Figure 10 and 
low dose rates in Figures 3-4, and vice versa. 

The baseline semi-empirical geomagnetic cutoff rigidity model will be based on lookup 
tables of a global grid of vertical cutoff rigidities calculated from the present to 2020. The 
vertical cutoff rigidities are calculated from numerical solutions of charged particle 
trajectories in the geomagnetic field, which is simulated using the International 
Geomagnetic Reference Field (IGRF) model [Barton, 1997], using the techniques 
advanced by Shea and Smart [1983] and Smart and Shea [2000, 2001], The IGRF model 
includes dipolar and non-dipolar contributions to the geomagnetic field. 

The global grid of vertical cutoff rigidities described above captures the cutoff rigidity 
during quiescent conditions. The simple Stonner rigidity relations [Stormer, 1930] for a 
dipolar field can be used to derive analytical scaling factors to adjust the tabulated 
vertical rigidities to allow for oblique angles of incidence of the GCR/SEP particles and 
the inclusion of geomagnetic storm effects [De Angelis et al., 2004; Wilson et al., 1991, 
2003a, b]. Accounting for oblique angles, by defining the cone of acceptance, is 
particularly important for the solid angle integration in HZETRN in calculating the 
particle flux and dose rates. 

The parameter that accounts for perturbations in the cutoff rigidities due to geomagnetic 
storms is a change in the horizontal component of the geomagnetic field at the magnetic 
equator [Kuhn et al., 1965; Wilson et al., 1991, 2003; De Angelis et al., 2004]. By 
applying a constant horizontal magnetic field to the particle trajectory equations describe 
by the Stonner theory, the vertical cutoff rigidity can be described by the following 
equation 


R = 14.9 cos 4 A 

m 



H.K { * 

M \ cos 6 X m 


\ 

1 

) 


(5) 


where R is the cutoff rigidity in units of GV. In the above equation, X m is the magnetic 
latitude, R e is the earth radius, M is magnetic dipole moment of the earth, and H st is the 
horizontal component of the storm-time magnetic field perturbation. Using the Dst-index 
to parameterize H st , Figure 1 1 shows the effect of geomagnetic field perturbation on the 

cutoff rigidity for recent solar-geomagnetic storms. The overall effect is a decrease in the 
cutoff rigidity, particularly at mid- to high-latitudes. Thus, particles that would normally 
be shielded (i.e., deflected) by the geomagnetic field during quiescent conditions have 
full access to the atmosphere during geomagnetic storm periods. 
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The semi-empirical geomagnetic cutoff model described above will be extended by 
driving it with real-time observations of the interplanetary magnetic field (IMF) and solar 
wind parameters. Real-time assessments of the change in the horizontal component of the 
geomagnetic field can be captured by dynamically including the effects of the ring 
current, which is included in the Tsyganenko [1989, 2002] T96 (and later models) 
empirical magnetic field model. 

In summary, the baseline geomagnetic cutoff rigidity model will consist of a tabulated 
global database of vertical cutoff rigidities computed from the present time period to 
2020 in one year increments using the techniques of Smart and Shea [2000, 2001], 
Analytic scaling relations based on Stormer theory will be used to define the angle of 
acceptance and to include dynamical perturbations to the horizontal component of the 
geomagnetic field [De Ange/is et a/., 2004; Wilson et al., 1991, 2003a, b]. Real-time 
perturbations to the horizontal component of the geomagnetic field will be modeled by 
driving the T96 magnetic field model [Tsyganenko, 1989, 2002] with measurements of 
IMF and solar wind pressure from the NASA/ACE satellite and Dst-index obtained from 
NOAA/SEC. 

A more sophisticated calculation of cutoff rigidity during magnetically disturbed periods 
can be provided by the use of a full numerical simulation using the coupled 
magnetosphere-ionosphere-thermosphere (CMIT) model [Wiltberger et al., 2004]. This 
model combines the Lyon-Fedder-Mobarry MHD simulation of magnetospheric 
dynamics [Lyon et al, 2004] with the Thermosphere-Ionosphere Nested Grid (TING) 
model for the upper atmosphere and ionosphere [ Wang et al., 1999]. It can be run using 
solar wind and IMF data such as is available in near-real-time from the NASA/ACE 
spacecraft. The Center for Integrated Space-weather Modeling (CISM) [ Luhmann et al, 
2004; Spence et al, 2004] is developing a comprehensive model for studying the 
interaction of solar energetic particles with the magnetosphere using a 3D Lorentz 
integration of SEP trajectories in electric and magnetic fields taken from the CMIT model 
[Kress et al, 2005, 2004; Hudson et al, 2004; Weygand and Raeder, 2005], It can be 
employed to obtain a detailed morphology of the cutoff rigidity using a dynamic 
magnetic field which results from geomagnetic activity that typically accompanies these 
events. We will further develop this capability in conjunction with the CMIT model and 
investigate its applicability to use with combined inputs from solar wind, IMF, and 
energetic particle data. These comprehensive physics-based geomagnetic cutoff 
simulations will also be used to benchmark the semi-empirical cutoff model and 
characterize its uncertainty. A common interface will be defined so that this 
comprehensive physics-based approach can be merged into the NAIRAS model as the 
MHD-based model matures and becomes computational feasible for real-time 
simulations. 

3.5 Meteorological Fields 

Both the AIR model and HZETRN are parameterized by atmospheric depth in units of 
g/cm . Real-time, global atmospheric depth is detennined by pressure and geopotential 
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height data obtained from the Met Office (METO) three-dimensional variational data 
assimilation (3-D-Var) system [ Lorenc et al., 2000], which is an update of the analysis 
method of Swinbank and O ’Neill [1994], Column abundance (or atmospheric depth in 
our units) is determined at each METO pressure surface by integrating atmospheric 
density over vertical height. Atmospheric depth (g/cm 2 ) is obtained at the commercial 
airline altitudes by linearly interpolating column densities at the METO pressure surfaces 
linearly in log pressure, using the geopotential height data at each pressure surface. 

Notice that the dose rates in Figures 1-4 are greater in the northern hemisphere polar 
region in January compared to July because the pressure level is lower in January, 
reducing the atmospheric shielding effect. 

4.0 Conclusions 

The Nowcast of Atmospheric Ionizing Radiation for Aviation Safety (NAIRAS) model is 
an analytical tool that will provide the first-ever, global, real-time, atmospheric ionizing 
radiation dosimetry package for archiving and assessing radiation exposure levels at 
commercial airline altitudes that have potentially harmful health outcomes. The planned 
web interface to NAIRAS will enable the radiation exposure levels of crewmembers to be 
monitored on a continual basis, as individuals or as groups. Flight path coordinates can be 
entered to track the radiation exposure for crewmembers as individuals. A “phantom” 
pilot feature will enable one to quickly evaluate the accumulated exposure levels for 
typical flight schedules on a representative set of flight paths (e.g., New York to London, 
Chicago to Hong Kong, etc.). NAIRAS data will provide the FAA and the commercial 
airline industry with valuable information for developing policies and procedures for 
modifying aircrew travel schedules so crewmembers do not exceed annual or career 
radiation exposure limits. SEP events are particularly worrisome since annual radiation 
exposure limits can be exceeded in one flight, and the potential for serious health 
outcomes is expected to be high with prenatal exposures especially risky. NAIRAS 
results will provide critical data for airline management decisions regarding flight-path 
alterations during SEP events, which must balance the cost incurred by rerouting against 
mitigating the potential health risks. NAIRAS results will also provide the input data for 
researchers and analysts in the public health sector (e.g., NIOSH) to improve 
understanding of the radiobiological effects and health outcomes of atmospheric ionizing 
radiation so refined policies can be established regarding radiation exposure limits. 
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Table 1. Current and Projected Maximum Allowable Exposure Limits (Wilson et al 1995) 




Maximum allowable exposure, mSv 


Exposure condition 

Present 
United States 
10CFR Part 20 
(1991) 

Proposed 
United States 
NUREG/BR-0 117 
(Cool and Peterson 1991) 

Proposed 
NCRP Rep. 
116 
(1993) 

Proposed 
ICRP Publ. 60 
(1991) 

Occupation: 

Annual 

Lifetime 

Pregnancy (total) 
Pregnancy (monthly) 

a 50 

[50 (Age - 18)] 
5 

50 

5 

50 

b 10 x Age 
0.5 

20 

c 2 

Public: 

Annual, many years 

d l 

1 

1 

1 

Annual, occasional 
Pregnancy (total) 
Pregnancy (monthly) 

5 

5 

5 

0.5 

c 2 


a Not to exceed 30 mSv in any quarter year. 
Recommended limit for new designs in 10 mSv/yr. 


c 

Abdomen surface for x-rays, 1 mSv in utero. 
d 5 mSv allowed with prior approval of NRC. 
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Figure 1: Neutron flux and ionization rate profiles for various latitudes. The blue lines represent 
solar maximum conditions. The green lines represent solar minimum conditions. The top row 
corresponds to January atmospheric conditions while the bottom row corresponds to July. 
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Figure 2: Neutron and ion dose-equivalent rates profiles at various latitudes. The blue lines represent 
solar maximum conditions. The green lines represent solar minimum conditions. The top row 
corresponds to January atmospheric conditions while the bottom row corresponds to July. 


20 


Dose Equivalent Rate for Solar Minimum at 12 km 
January mSv/IOOO hr July 



Figure 3: Northern-hemisphere dose-equivalent rates at 12 km for solar minimum conditions. The 
atmospheric condition for the left figure is January and the right figure is July. 
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Dose Equivalent Rate for Solar Maximum at 12 km 
January mSv/IOOO hr July 



Figure 4: Northern-hemisphere dose-equivalent rates at 12 km for solar maximum conditions. The 
atmospheric condition for the left figure is January and the right figure is July. Note the difference 
in scale compared to Figure 3. 
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Figure 5: Differential flux computed by HZETRN at various atmospheric depths for solar minimum 
conditions. The figures show neutron flux and selected ion flux. 
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Figure 6: Total dose-equivalent rate computed by HZETRN at various atmospheric depths for solar 
minimum conditions. 
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Figure 7: Event-integrated differential proton and neutron flux computed by HZETRN at various 
atmospheric depths for the September 1989 SEP event. 
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Figure 8: Event-integrated total dose-equivalent rate computed by HZETRN at various atmospheric 
depths for the September 1989 SEP event. 


26 




Baddwar and QNelll GCK Model 



Figure 9: GCR spectra for selected ions. The solid-red line is the local interstellar spectrum (LIS). 
The dashed-green line corresponds to solar minimum conditions while the dotted-blue line 
corresponds to solar maximum conditions. 
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Figure 10: Northern-hemisphere vertical geomagnetic cutoff rigidities computed from particle 
trajectories using the IGRF model. 
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Vertical Geomagnetic Cutoff Rigidity- 



Figure 11: Zonal-averaged vertical geomagnetic cutoff rigidity. The solid-red line corresponds to the 
quiescent cutoff rigidities computed from particle trajectories and the IGRF model. The blue lines 
correspond to cutoff rigidities computed from (5) using the Dst-index to parameterize the horizontal 
magnetic field perturbation for recent solar-geomagnetic storms. 


29 


